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Pb 2 V3 0g is a quasi one-dimensional compound with antiferromagnetic alternating chains. 
We measured the magnetization and the specific heat in a single crystal of PD2V3O9 with 
applied magnetic fields up to 14 T and established the magnetic phase diagram. The phase 
boundary in the high magnetic field region is in a good agreement with that predicted by 
the mean field theory in the classical magnetic ordering of the transverse magnetization, 
while that in the low magnetic field region should be discussed in the viewpoint of the Bosc- 
Einstcin condensation of the field induced magnons. Thus, the experimentally determined 
phase diagram shows the asymmetry which should be originated in the mixing contribution 
of excited states as well as in the interchain interactions. 

KEYWORDS: Pb 2 V30g, Bose-Einstein condensation, magnons, alternating chain, field in- 
duced magnetic ordering, magnetic phase diagram 

1. Introduction 

Low dimensionality in quantum spin systems has been a matter of particular interest 
over some decades since it could lead to novel exotic phenomena. Since the ground state in 
the low-dimensional system has instabilities with quantum fluctuations, it can be controlled 
by external fields or impurity effects. A weakly coupled spin gap system is a candidate for 
exhibiting novel quantum phenomena due to its low dimensionality. In the case of dimer 
systems without any interdimer interactions, a dimer has the ground state of \S,S Z ) = |0,0) 
with excited triplet states of \S, S z ) = |1, 1), |1, 0), |1, — 1). The field dependence of the energy 
levels is shown in Fig. 1 as solid lines. The external magnetic field splits the degenerated triplet 
states and reduces the energy gap between the |0, 0) and |1, 1} states. Under the magnetic field, 
the ground state changes from the |0, 0) to 1 1, 1) state and spins induced by the field are aligned 
along the field. In the presence of the additional antiferromagnetic interdimer interactions, 
the field-induced magnetic moments can exhibit the magnetic ordering of their transverse 
spin components. The antiferromagnetic long range ordering (AFLRO) in this manner was 
predicted by the mean field theory. ^ The AFLRO state was approximated by the ground 
state of mixed |0, 0) and |1, 1) states. The energy levels of the mixed states derived based on 
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the mean-field magnetic ordering model of the transverse magnetization are shown in Fig. 1 
by dashed lines. In addition, the symmetric bell-shape phase boundary can be derived. 1 ' The 
lower critical field H c i, the upper critical field H C 2 and the maximum transition temperature 
to AFLRO phase T max can be calculated from the intradimer and interdimer interactions. 
At low temperatures, the phase boundary becomes almost perpendicular to H axis. In fact, 
TlCuCl 3 , 2 ' 3 ) KCuCl 3 , 4 ) Ba 3 Cr 2 8 , 5 ) Sr 3 Cr 2 8 6 ) and Pb 2 V 3 9 7 - 10 ) were found to show the 
field-induced magnetic orderings. The induced |1,1) excitations would have characteristics 
of bosons as excited magnons. The AFLRO phase can be interpreted in terms of a Bose- 
Einstein condensation (BEC) of magnons. The increase of the magnetic susceptivility below 
the transition temperature is explained by the increase of condensed magnon density. The 
magnon BEC picture also gives a power-law like critical behavior described as 

T c « [H cl -H cl {T = Q)] 1 't (1) 

where T c is the transition temperature to AFLRO phase and <j> is the critical exponent. 
Most experimental data supported the magnon BEC picture. The critical exponent 4> was 
experimentally reported to be between 1.5 and 2.3, while the value was theoretically predicted 
as 1.5. 2 ) Furthermore, from the view of the phase diagram topology, asymmetric H — T phase 
diagrams was reported in S^C^Os 6 * 1 and Pb 2 V 3 0g, 8 ) while symmetric one was expected 
because of the particle - hole symmetry even in the case with the BEC model. 

Pb 2 V 3 0g was first synthesized by the electrochemical reduction of molten PbV 2 06 11 ^ and 
found to have antiferromagnetic alternating chains. The arrangement of the chains in Pb 2 V 3 0g 
is shown in Fig. 2(a) and the alternating chain is illustrated in Fig. 2(b). A magnetic V 4+ ions 
with coordinating six oxygen atoms form V 4+ 06 octahedrons with corner-sharing of oxygen 
atoms, forming the spin chain with S = 1/2 along the [101] direction. The octahedrons are 
aligned on a tilt and the distances between adjacent V 4+ ions change alternatively. Along 
the [101] direction, V 5+ 04 tetrahedrons stand between V 4+ 06 octahedrons. The calculation 
by the quantitative spin dimer analysis based on tight binding calculations suggested that 
the V 4+ -0-0-V 4+ super-superexchange interactions along the [101] direction were dominant 
for the antiferromagnetic alternating chain model. 12 ) From magnetization measurements, the 
existence of a spin gap and the phase transition to the AFLRO phase under magnetic fields 
have been clarified. 7 ) Lower and upper critical fields have been determined as H c \ = 4 T and 
H C 2 = 38 T, respectively. From the power low behavior of phase boundary, the critical expo- 
nent has been determined as 4> = 1.9, 7 ) or 2.O. 9,10 ) In this paper, we report the magnetic phase 
diagram precisely determined by using a single crystal in order to elucidate the asymmetric 
H — T phase diagram. 
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2. Experiments 

The polycrystalline sample of PD2V3O9 was prepared by solid state reaction from PbO, 
V2O5 and V2O3. The powders were mixed and then sintered at 570 C° for 2 days with in- 
termediate grindings. Single crystals of PD2V3O9 were grown by self-flux and floating zone 
methods. The crystal growth by self- flux method was made by cooling the sintered sample 
slowly from 750 C° in an evacuated quartz tube. Glossy crystals were obtained with a max- 
imum size of 0.2 x 0.5 x 0.2 mm 3 . By floating zone method, 8 g of powdered sample was 
molded into a rod under pressure of 400 kgf/cm 2 . The rod was sintered at 500 C° for two 
days in vacuum and then was set to the sample holder under Ar atmosphere of 1 atm. The 
crystal growth was made by melting the sample rod partially with a crystal growth speed of 
0.5 mm/h. The size of obtained glossy crystal is about 2x2x1 mm 3 . A Debye-Scherrer 
X-ray diffraction measurement was performed to check the purity, and the single crystal X- 
ray diffraction measurement was done by transmission Laue method to determine the crystal 
plane, and axes. The magnetic susceptibility was measured by a superconducting quantum 
interference device magnetrometer. The magnetization curve was measured up to 55 T with 
a pulse magnet at Ultra High Magnetic Field Laboratory, Institute for Solid State Physics, 
University of Tokyo. The specific heat was measured by relaxation method under magnetic 
field up to 14 T. 

3. Results 

Viewgraphs of single crystals grown by floating zone method are shown in Fig. 3(a). Figure 
4 shows the Debye-Scherrer X-ray diffraction patterns. The crashed powders of single crystals 
were confirmed to be in a single phase. In addition, each single crystal was found to have a 
single domain by transmission Laue X-ray photographs. The crystal axes of the sample used 
in measurements of physical properties were determined as shown in Fig. 3(b). 

Figure 5(a) shows the temperature dependence of the magnetic susceptibility under the 
magnetic field of 1 T. The susceptibility does not depend on the field direction so much. This 
fact suggests that the magnetic anisotropy is not so large that the orbital contribution would 
be small. To estimate spin interactions, the susceptibility was fitted to the isolated spin dimer 
model 13 ) as well as to the antiferromagnetic alternating chain model. 14 ) The fitting functions 
are described by 

C 

Xo = Xspin(#, Jo, Jl) + ^ + Xo, (2) 

where Jo, Ji,g, C, xo are the intradimer interaction, the interdimer interaction, g-factor, the 
Curie constant and the temperature independent term, respectively. The values of estimated 
parameters are shown in Table I. In the case of the isolated dimer, of course, J\ must be zero. 
The alternating chain model reproduces the experimental data much better than the isolated 
dimer model does, indicating the presence of the interdimer interaction. Here, the Jq and J\ 



3/14 



J. Phys. Soc. Jpn. Full Paper 

were determined as 31.5 and 23.0 K, respectively. 

The high- field magnetization curve measured at 1.3 K is shown in Fig. 6. The whole be- 
havior is almost consistent with that observed in the previous report using the polycrystalline 
sample. 8 ^ The critical fields in 1.3 K were determined as H c \ = 4.7 T and H C 2 = 38.4 T from 
anomalies in the field derivative of the magnetization curve. 

The specific heat was measured to determine the phase diagram precisely. Figure 7 shows 
the specific heat at temperatures between 0.5 and 6 K with the field applied along or perpen- 
dicular to the b* axis. Below 4.5 K, a A-type anomaly due to the second order phase transition 
was observed, which corresponds to the field induced magnetic ordering. With increasing field, 
the peak shifts to the higher temperature and the intensity is enhanced. The shift of peak by 
the magnetic field can be explained by the description of BEC of magnons. The critical tem- 
perature T c was determined from the second temperature-derivative of {C(H) — C(H = 0)}/T, 
where C(H) is the specific heat under the field H. Rough estimation of the entropy changes 
accompanied by the phase transition was given by integrating C{H)/T after extracting a 
phonon contribution of Cphonon ~ 2.94 T 3 mJ/(K • mol) from C(H), as shown by solid line in 
Fig. 7. The phonon estimation was made with an assumption that the entropy change in the 
alternationg chain model agrees with that in the isolated dimer model at high temperatures. 
The Jo in Table I was used for the entropy estimation of the isolated dimers. Figure 8 shows 
the estimated entropy with the field applied along the b* axis. The entropy estimation with 
the field perpendicular to b* axis provides a quite similar result. The entropy is increased 
with increasing the field and suppressed below T c . The enhancement by the field results from 
the reduction of the spin gap and the suppression below T c indicates the magnetic ordering 
under the magnetic field. In the field of 6 T, the entropy is estimated as 0.12 J/(K • mol), 
which corresponds to the 2.1 % of the maximum value i?ln2 = 5.76 J/(K mol) expected for 
the spin system of S = 1/2. Therefore, the critical magnon density in 6 T [T c = 2.1) can be 
described as n c = 0.021 and becomes about 10 times larger than that in TlCuCls 3 * 1 at the same 
temerature. This indicates that more magnons are induced in Pb2V30g than in TICUCI3. 

The phase diagram was determined as shown in Fig. 9. Whole behavior of the phase 
boundary does not depend on the field direction. The anisotropic interactions could yield 
the field-dependence of the phase diagram, and in Pb2VsOg, e.g., the Dzyaloshinskii-Moriya 
interaction could work due to the absence of the inversion symmetry at the center of V 4+ 
pairs. In the case of Pb2V30g, however, the weak field-direction dependencies in x, C and the 
phase diagram of Pb2VsOg indicate that the anisotropic interactions are small. 

4. Discussion 

Here we discuss the shape of the phase diagram and its asymmetric behavior in this section. 
The dashed line in Fig. 9 expresses the phase boundary of the classical transverse magnetic 
ordered phase, which can be produced by the mean field theory with the interaction parameters 
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estimated from the x ~ T curve. The spin Hamiltonian of one-dimensional alternating chain 
is described by 

H = H + H' 

H = Jo^2 Sji • s i2 - gHB~R ■ ^2(sn + 8 i2 ) . . 

i i ^ ' 

H' = J\ ^ Sj2 • Sj+l l, 

i 

where i represents the notation of the dimer and il and i2 denote the two spins on each dimer. 
The site i and i+1 locate on the nearest neighbor along the spin chain. The arrangement of 
the interaction parameters and the notations are shown in Fig. 2(b). The direction of the spin 
chain is not important in the following discussion. According to the mean field theory, Tachiki 
and Yamada confined themselves to subspace of the |1, 1) and 1 0, 0) states in the vicinity of 
H c \}^ The model Hamiltonian in eq. (3) is rewritten into the following low energy effective 
Hamiltonian as, 

H = H + H' 

= T ( J ° + Y ~ 911bH ) 1 ~ \ ( J ° + T ~ 9 ^ bH ) aiz ( 4 ) 

+ ^ 2(cTi x cr i+ i x + ai y a i+ i y ) + Ui z Gj z ]. 

i 

From eq. (4), H c i, and T max can be estimated as 

H cl = -L(j -±f)= 14.2 T 
»\ 2 ) 

H c2 = — (J + Ji) = 38.8 T (5) 

9VB 

T max = = 5.8 K. 
Ak B 

The estimated value of H c2 is quite close to the experimental one. The magnetic field at which 
T c becomes the highest is estimated to be H = (H c \ + H c2 )/2 = 26.5 T, which is also quite 
close to the experimental value as shown in Fig. 9. In contrast, the estimated H c \ deviates 
largely from the experimental value. 

Figure 10 shows the phase boundary focused on the lower critical field H c ±. The power law 
behavior of eq. (1) was predicted in the magnon BEC model. The dashed and solid lines in 
Fig. 8 are the phase boundaries produced by the fitting to eq. (1). The dashed line is made by 
fitting the experimental points under 6 K with cp confined to the value of 1.5 and the solid line is 
made by using 4> as a fitting parameter, which is determined as 2.1. While the critical exponent 
4> is expected to be 1.5 in the magnon BEC model, the experimental data are better explained 
by cj) = 2.1 with H c \ = 3.4 T. In early papers, the critical exponent exhibits (j) = 1-9, 7 - ) or 
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2.0 9 ' 10 ) and the possibility of <p convergence to the universal value 4> = 1.5 was suggested in 
the low-temperature limit, which was first observed in the experimental phase boundary of 
TlCuCbj. The BEC theory by using realistic dispersion relations of TlCuCbj explained the 
phase boundary below 5 K and <j) = 1.5 was obtained from fitting in the temperature range 
below 2.4 K. 3 ) Also in the case of Pb2V30g, the valid description of the magnon dispersion 
should be necessary for the universality of AFLRO phase. 

To evaluate the magnon-magnon interactions more correctly, we took into consideration 
the additional effects of the |1,0) and |1,— 1) states. In the subspace of the 1 0, 0} and |1,1) 
states, a point inversion symmetry h e R — > —h e g,a — > — a in Hamiltonian eq. (4) ensures the 
symmetry of the phase diagram for a zeropoint of the effective field h e s = h — Jq — Ji/4. 
Second-order processes between the |0,0), |1,1) states and the excited 1 1 , 0) , |1,— 1) states 
yield the perturbation terms and they break the inversion symmetry in the spin Hamiltonian. 
The mixing effect with the excited states was discussed in Refs. 15 and 16. In these references, 
a bond-operator technique was used to describe both the spin-liquid and the ordered states in 
the field and pressure induced quantum phase transition in TICUCI3. The values of H c \ and 
H C 2 were described by intradimer and interdimer interactions, and magnetization curves and 
the magnon dispersion relations observed in TICUCI3 were accounted for by the bond-operator 
technique. In this theory, the operators which correspond to creations and annihilations of 
the |0, 0), |1, 1), |1, 0} and |1, —1) states represent spin degrees of freedom on each dimer. The 
interdimer interactions yield the mixed ordered states which are approximated by a linear 
combination of the |0, 0), |1, 1} and |1, —1) states. The |1,0) state does not contribute its term 
to the mixed states. In consideration of the mixing contribution of the |1, — 1) state with the 
|0, 0) and |1, 1) states, H c \ of Pb2V30g is obtained as 

H cl = Jo- = H.7 T (6) 

and H C 2 becomes exactly the same as that in eq. (5). The in eq. (5) agrees with the 
experimental H C 2. The H c \ in eq. (6) is smaller than that in eq. (5) and becomes closer to 
the experimental H c ±. The asymmetric behaviors especially in the vicinity of H C 2 and H c \ 
are due to the change in mixing contribution of the |1,— 1) state. Around H C 2, the mixing 
contribution is ignorable and therefore the magnon-magnon interactions in the high field 
region can be described by quasi one-dimensional mean field. In the vicinity of H c i, the 
finite mixing contribution reduces the minimum of the magnon dispersion and makes H c \ 
smaller than the value in eq. (5). In addition, magnon-magnon interactions are well described 
by interacting dilute bosons and therefore the phase boundary shows power-law behavior. 
The alternating chain model with the field-induced magnetic ordering explains the physical 
properties of Pb2V30g very well. Details of AFLRO state in Pb2V30g, especially the critical 
exponent, may be explained with the modified alternating chain model by additional small 
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interchain interactions. The experimental H c \ is reproduced well if the value of J\ is regarded 
as 30.5 K. Since the assumed J\ is not so different from that estimated from X~T curve, the 
presence of interdimer interactions may increase J\ and also affect the critical exponent. 

5. Conclusion 

In summary, the magnetic phase diagram of PD2V3O9 was determined precisely from the 
high field magnetization and the specific heat by using single crystal. The experimental phase 
boundary near H& coincides with that based on the mean field picture, while that in the 
vicinity of H c \ shows the power-law behavior predicted in the magnon BEC picture. These 
results indicate that the mixing contribution of excited states with low-lying states shifts the 
phase boundary to the lower field especially in the case of H& . 
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II, 1> 



Fig. 1. (Color online) Energy levels of the spin states of \S, S z ) = |0, 0), |1, 1), |1, 0) and |1, —1) in the 
dimers without any interdimer interactions (solid lines) and the mixed \S, S z ) = |0, 0) and |1,1) 
states in interacting dimers derived by the mean field theory (dashed lines). 1 ) 




Fig. 2. (Color online) (a) The arrangement of V 4+ Og octahedrons and V 5+ C>4 tetrahedrons in the ac 
plane, (b) The illustration of the alternating chain with the notations above and the interaction 
parameters below the chain, both of which are used in eq. (3) (see text). 
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Fig. 3. (Color online) (a) Viewgraphs of single crystals obtained by floating zone method, (b) Crystal 
axes of the single crystal determined by transmission Laue method. 
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Fig. 4. (Color online) Debye-Scherrer X-ray diffraction patterns for the single and powdered crystals 
with calculated powder pattern for comparison. 
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Pb 2 V 3 9 H = 1 T 
A H±a, HUb* 
□ Hlla 
• Hllb* 




Fig. 5. (Color online) (a) Temperature dependence of magnetic susceptibility under the applied field 
of 1 T parallel and perpendicular to a, b* axes, (b) Temperature dependence of magnetic suscep- 
tibility with H || b* (closed circles). Dashed and solid lines show the fitted lines with the isolated 
dimer 13 ) and the alternating chain 14 ) models, respectively. 
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Fig. 7. (Color online) Temperature dependences of specific heats with the field (a) H \\ b* and (b) 
H _L b* , with offset of H x 0.1 J/(K • mol). Solid line shows the phonon contribution. 
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Fig. 8. (Color online) Temperature dependence of magnetic entropy with the field H \\ b* . Solid arrow 
indicates T c at each field. Solid line in inset shows the magnetic entropy of the isolated dimer with 
the spin gap A = 26.7 K under zero field. 
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x re f- 8) phase boundary 

• T c (EC,Hllb*) deduced from 
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■ T c (HC,HXb*,HXb*) approximation 
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Fig. 9. (Color online) The H — T magnetic phase diagram of Pb2V 3 0g. Solid line is a guide for eyes. 
Dashed line corresponds to deduced phase boundaries from classical mean field approximation 1 ' 
based on the parameters in Table I. 
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Fig. 10. (Color online) Magnetic phase diagram focused on at low temperatures and in the low field 
region. Dashed and solid lines are obtained by power law as T = A(H — H c x(T = O)) 1 ^ with 
critical exponent tj> = 1.50 and <j) = 2.09, respectively. 
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Table I. Parameters obtained from fitting the magnetic susceptibility to the isolated dimer 13 ) and 
the alternating chain 14 ) models. 





Isolated dimer 


Alternating chain 


Jo (K) 


26.7 ±0.3 


31.5 ±0.1 


Ji (K) 




23.0 ±0.9 


g 


1.75 ±0.02 


2.09 ±0.01 


C (emu • K/mol) 


(6.78 ±0.47) x 10~ 3 


(5.86 ±0.13) x 10~ 3 


Xo (emu/mol) 


(6.04 ±0.61) x 10~ 4 


(-1.08 ±0.18) x 10- 4 
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